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Abstract 

We provide nonperturbative fragmentation functions for B mesons, both at lead- 
ing and next-to-leading order in the MS factorization scheme with five massless 
quark flavors. They are determined by fitting the fractional energy distribution of 
B mesons inclusively produced in e + e~ annihilation at CERN LEP1. Theoretical 
predictions for the inclusive production of B mesons with high transverse momenta 
in pp scattering obtained with these fragmentation functions nicely agree, both in 
shape and normalization, with data recently taken at the Fermilab Tevatron. 
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1 Introduction 



The study of 5-quark production in high-energy hadronic interactions offers the opportu- 
nity to test perturbative quantum chromodynamics (QCD) Hadron-collider experi- 
ments usually consider the cross section integrated over a fixed range in rapidity r\ and 
over all values of transverse momentum p? above a variable threshold p™ m . First mea- 
surements of such cross sections were performed by the UAl collaboration at the CERN 
pp collider with center-of-mass (CM) energy yfs = 630 GeV ||. More recent experimental 
results at y/s = 1.8 TeV were presented by the CDF || and DO collaborations at the 
Fermilab Tevatron. On the theoretical side, such cross sections were calculated up to 
next-to-leading order (NLO) in the strong coupling constant a s |1]J5|,|6|]. The shapes of the 
theoretical curves agree well with the data of all three experiments, UAl, CDF, and DO. 
However, independent of the beam energy, the absolute normalizations of the experimen- 
tal cross sections exceed, by about a factor of two, the respective predictions obtained 
with the conventional scale /i = \jp\ + ml and a typical 6-quark mass of m& = 4.75 GeV. 
The experimental cross sections could only be reproduced by the theoretical predictions if 
/i and m b were reduced to /i = Jp\ + ml/2 and m b = 4.5 GeV and parton density func- 
tions (PDF's) with particularly large values of the asymptotic scale parameter A were 
chosen. 

In the experimental studies, the 6-quark production cross section is not actually mea- 
sured directly as a function p™ m . The original measurements refer to the production of B 
mesons, which decay either semileptonically, or exclusively or inclusively into J/ip mesons. 
The cross sections for the production of bare h quarks were then obtained by correcting 
for the fragmentation of b quarks into B mesons with the help of various Monte Carlo 
(MC) models. Since this is a model-dependent procedure, it remains unclear whether the 
disagreement between the experimental data and the NLO predictions is actually real. In 
order to extract the 6-quark production cross section, one needs an independent measure- 
ment of the fragmentation of b quarks into B mesons. In fact, two years ago the CDF 
collaboration published data on their first measurement of the 5-meson differential cross 
section da/dpT for the exclusive decays B + — > J/tpK + and B° — > J/tpK*° based on the 
1992-1993 run (run 1A) 0. Similarly to the case of 6-quark production, the measured 
cross section was found to exceed the NLO prediction by approximately a factor of two, 
while there was good agreement in the shape of the pt distribution. Agreement in the 
normalization could only be achieved by choosing extreme values for the input parameters 
of the NLO calculation, i.e., by reducing /i and m b and by increasing A 0. 

In Ref. ||, the CDF collaboration extended their analysis by incorporating the 
data taken during the 1993-1996 run (run IB), which yielded an integrated luminosity of 
54.4 pb _1 to be compared with 19.3 pb _1 collected during run 1A ||, and presented the 
cross section da/dpT for the inclusive production of B + and B° mesons with p? > & GeV 
in the central rapidity region < 1. Again, the NLO prediction with input parameters 

similar to those used for the integrated 6-quark cross section in Refs. [l],§,[| (/i = \jp\ + ml 
and m& = 4.75 GeV) was found to agree with the data in the shape, while its normalization 
came out significantly too low, by a factor of 2.1 ± 0.4. Here, it was assumed that the 
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fragmentations of b quarks into B mesons can be described by a Peterson fragmentation 
function (FF) || with e = 0.006. This value for the e parameter was extracted more than 
ten years ago from a global analysis of data on 5-meson production in e + e~ annihilation at 
PEP and PETRA |L0|, based on MC models which were then up-to-date. The result of this 
comparison is to be taken with a grain of salt, since the underlying description of b — ► B 
fragmentation is ad hoc and not backed up by model-independent data. It is the purpose 
of this work to improve on this situation. The fragmentation of b quarks into B mesons 



has been measured by the OPAL collaboration at LEP1 The produced B + and 

B° mesons were identified via their semileptonic decays containing a fully reconstructed 
charmed meson. This resulted in the measurement of the distribution of the B mesons in 
the scaling variable x = 2E(B)/y/s, where E(B) is the measured energy of the B + /B° 
candidate and yfs = Mz is the LEP1 CM energy. Earlier measurements of the b — > B FF 



were reported by the L3 collaboration at LEP1 [fL^j . In the following, we shall base our 
analysis on the OPAL data, which have higher statistics and contain more x bins than 
the L3 data. 

At LEP1, B mesons were dominantly produced by Z — > bb decays with subsequent 
fragmentation of the b quarks and antiquarks into B mesons, which decay weakly. In the 
reaction e + e~ — > bb — > B + X at the Z-boson resonance, the b quarks and antiquarks 
typically have large momenta. A large-momentum b quark essentially behaves like a 
massless particle, radiating a large amount of its energy in the form of hard, collinear 
gluons. This leads to the presence of logarithms of the form a s ln(M|/m^) originating 
from collinear singularities in a scheme, where m& is taken to be finite. These terms 
appear in all orders of perturbation theory and must be resummed. This can be done by 
absorbing the m^-dependent logarithms into the FF of the b quark at some factorization 
scale of order Mz- Alternatively, one can start with = and factorize the collinear 
final-state singularities into the FF's according to the MS scheme, as is usually done in 
connection with the fragmentation of light quarks into light mesons. This is the so-called 
massless scheme [O], in which m& is neglected, except in the initial conditions for the FF's. 



This scheme was used for NLO calculations of charm and bottom production in e + e [14 



pp [T3], 7J9 [IB, 17 1, and 77 Wjfy collisions, with the additional feature that the massless c 
and b quarks were transformed into their massive counterparts by the use of perturbative 
FF's ||14|| . These perturbative FF's enter as a theoretical input at a low initial scale /j,q 
of order m c or mj, respectively, and are subject to evolution to higher scales /i with the 
usual Altarelli-Parisi (AP) equations fl9|]. Following Ref. p0 |, this theory was extended 
by including nonperturbative FF's, which describe the transition from heavy quarks to 
heavy mesons jl7| , |2ll , p2| . 

In this work, we describe the fragmentation of massless b quarks into B mesons by 
a one-step process characterized entirely in terms of a nonperturbative FF, as is usually 
done for the fragmentation of u, d, and s quarks into light mesons. We assume simple 
parametrizations of the 6-quark FF at the starting scale. We determine the parameters 
appearing therein through fits to the OPAL data |TT[ at lowest order (LO) and NLO. 



These 6-quark FF's are then used to predict the differential cross section da/dpx of B- 
meson production in pp scattering at y/s = 1.8 TeV, which can be directly compared with 
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recent data from the CDF collaboration ||. 

This paper is organized as follows. In Sect. 2, we recall the theoretical framework for 
the extraction of FF's from e + e~ data, which was previously used for c-quark fragmenta- 
tion into D*^ mesons | 23fl , and present our results for the 6-quark FF's at LO and NLO in 
the MS factorization scheme with five massless flavors. We assume three different forms 
for the FF's at the starting scale, which enables us to assess the resulting theoretical 
uncertainty in other kinds of high-energy processes, such as pp scattering. In Sec. 3, we 
apply the nonperturbative FF's thus obtained to predict the cross section of 5-meson 
production in pp collisions at the Tevatron and compare the result with recent data from 
CDF IHI. Our conclusions are summarized in Sec. 4. 



2 .B-meson production in e + e collisions 

Our procedure to construct LO and NLO sets of FF's for B mesons is very similar to the 
case of D*^ mesons treated in Refs. |]22|j23[| . Here we only give those details which differ 



from Refs. 22,23 



The OPAL data on the inclusive production of B + and B° mesons in e + e~ annihilation 
at the Z-boson resonance serve as our experimental input |]TT|. These data are presented 
as differential distributions in x = 2E(B)/ \fs, where E(B) is the measured energy of 
the B + or B° candidate. This function peaks at fairly large x. For the fitting procedure 
we use the experimental x bins, with width Ax = 0.08, in the interval 0.28 < x < 1 
and integrate the theoretical functions over Ax, which is equivalent to the experimental 
binning procedure. There is a total of nine data points. 

When we talk about the b —> B FF, we have in mind the four fragmentation processes 
b — > B + , b — > B°, b — > B~ , and b — > B°. In Ref. ||24}| , the respective branching fractions 



are all assumed to be equal. If we neglect the influence of the electroweak interactions, 
this follows from the u d flavor symmetry and the charge-conjugation invariance of 
QCD. We thus make the stronger assumption that the FF's of these four processes all 
coincide. We take the starting scales for the FF's of the gluon and the u, d, s, c, and b 
quarks and antiquarks into B mesons to be fiQ = 2m;,, with mj = 5 GeV. The FF's of the 
gluon and the first four quark flavors are assumed to be zero at the starting scale. These 
FF's are generated through the \i evolution. For the parametrization of the 6-quark FF 
at the starting scale /^o, we employ three different forms. The first one is usually adopted 
for the FF's of light hadrons, namely 

D b (x, f i ) = Nx a (l-x) fS . (1) 



This form has been used in Ref. [z(J to describe the nonperturbative effects of 6-quark 
fragmentation, in addition to a perturbative contribution. The standard (S) parametriza- 
tion ([Tp depends on three free parameters, N, a, and (3, which are determined by fits to 
the OPAL data [|ll[] after evolution to the factorization scale Mf = Mz- As our second 
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parameterization, we use the Peterson (P) distribution ||, 

^(*./»)= V- (1 ^+V (2) 

This choice is particularly suited to describe a FF that peaks at large x. It has been 
frequently used in connection with the fragmentation of heavy quarks, such as c or b 
quarks, into their mesons. It depends only on two parameters, iV and e. 

The third parametrization is theoretically motivated. There exists a particular class of 
FF's which are calculable in perturbative QCD, namely those of gluons and heavy quarks 
into heavy- heavy bound states, such as cc, bb [53], and cb mesons These perturbative 
FF's can also be applied to describe the fragmentation of b quarks into bound states of 
b and light quarks, in the sense of a model assumption rather than a formula derived in 
perturbative QCD. The formula for the b — > B c transition was derived by Braaten (B) et 
al. p6[ and reads 

D b (x,n ) = N— K - J — 6 - 18(1 - 2r)x + (21 - 74r + 68r 2 )x 2 

[1 — (1 — r)x\ b 1 

- 2(1 - r)(6 - 19r + 18r 2 )x 3 + 3(1 - r) 2 (l - 2r + 2r 2 )x 4 ] , (3) 

where r = m c /(mb + m c ) and N is given in terms of a s , m c , and the I? c -meson wave 
function at the origin. Similar formulas also exist for b — > B*,B** p6[ . Naively applying 



this formula for r to the fragmentation process b — > B would yield a rather small number, 
which is not well determined. Thus, our philosophy is to treat N and r as free parameters if 
one of the quarks in the bound state is light. In Ref. [g(J, the branching fraction of c — > B c 
was found to be two orders of magnitude smaller than the one of b — ► B c . Extrapolating 
to the case of B mesons, it hence follows that our assumption D q (x,fio) = 0, where q 
denotes a light quark, should be well founded even if q is the light constituent of the B 
meson. 

We calculate the cross section (l/<r tot )da/dx for e + e~ — > 7, Z — > B + / B° + X to LO 



and NLO in the MS scheme with five massless quark flavors as described in Ref. |27| , 
where all relevant formulas and references may be found. In particular, we choose the 
renormalization and factorization scales to be \x = Mf = ^/s. As for the asymptotic 
scale parameter appropriate for five active quark flavors, we adopt the LO (NLO) value 



Ag) = 108 MeV (227 MeV) from Ref. 0. As in Ref. J22J, we solve the AP equations in 



x space by iteration of the corresponding integral equations. In the Appendix of Ref. [22 



the timelike splitting functions are listed in a convenient form, i.e., with the coefficients 



of the delta functions and plus distributions explicitly displayed. As in Ref. |p2[ , we take 
the 6-quark mass to be = 5 GeV. Since only enters via the definition of the starting 
scale /j,Q, its precise value is immaterial for our fit. 



The OPAL data are presented in Fig. 3 of Ref. |L1[ as the distribution dN/dx nor- 
malized to the bin width Ax = 0.08. In order to convert these data to the inclusive 
cross section (l/<j tot )da/dx, we need to multiply them by the overall factor 2Rbf(b — ► 
B)/Ax = 2.198, where Rb = T(Z — ► bb)/T(Z — > hadrons), f(b —>■ B) is the measured 
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b — > B branching fraction, and the factor of two accounts for the fact that our cross- 
section formula ]27| includes the fragmentation of both b and b. Following Ref . [ 24 ] , we 



identify f(b — ► B) = f(b — > B + ) = f(b — > B°). For consistency, we adopt the OPAL 
results R b = 0.2171 ± 0.0021 ± 0.0021 H and f(b -> B) = 0.405 ± 0.035 ± 0.045 ||, 
where the first (second) error is statistical (systematic). 

The values for the input parameters in Eqs. ([[]), (@), and (0) which result from our 
LO and NLO fits to the OPAL data are summarized in Table [l]. In the following, we 
refer to these FF's as sets LO S, NLO S, LO P, NLO P, LO B, and NLO B, respectively. 
The corresponding \ 2 values per degree of freedom (xdf) are listed in the last column of 
Table [I]; there is a total of nine degrees of freedom. Except for the sets of type S, the 
Xdf values for the NLO fits are slightly lower than those for the LO fits. The Peterson 
ansatz (0) yields the best fits. This is surprising, since it has only two free parameters, 
one less than the the standard form ([J). The sets of type B have the largest Xdf vam es. 
Since the 6-quark FF is peaked at x ^> 0.5, we have a 3> ft in the case of sets LO S and 
NLO S. The e parameters of sets LO P and NLO P are larger than the standard value 
e = 0.006 []n| usually quoted in the literature. It is important to note that the values 
of e obtained in the various analyses depend on the underlying theory for the description 
of the fragmentation process b — > B, in particular, on the choice of the starting scale 
/j,q, on whether the analysis is done in LO or NLO (as may be seen from Table [l|), and 
on how the final-state collinear singularities are factorized in NLO. We emphasize that 
our results for e in Table [1] refer to the pure MS factorization scheme with five massless 
flavors and //o = Ivn^ = 10 GeV. If we were to interpret the values for r in Table [l] with 
the formula r = m q /(mb + Tn q ), which is naively adapted from the analogous definition 



for cb bound states p6fl , then we would find m q = 688 MeV and 924 MeV at LO and 



NLO, respectively. These values are a factor of 2-3 larger than the generally assumed 
constituent- quark masses of the u and d quarks. This just illustrates the model character 
of using ansatz (|||) in connection with heavy-light bound states. 

In Figs. |l|(a)-|l](c), we compare the OPAL data [T]J with the theoretical results eval- 



uated with sets S, P, and B, respectively. Except at low x, the LO and NLO results are 
very similar. At low x, we observe significant differences between LO and NLO. In this 
region, the perturbative treatment ceases to be valid. Here, the massless approximation 
also looses its validity. Since B mesons have mass, m(B) = 5.28 GeV, they can only be 
produced for x > x m j n = 2m(B)/Mz = 0.12. The LO result has a minimum in the vicinity 
of £ m m and strongly increases as x — > 0. Therefore, our results should only be considered 
meaningful for x^x cut with x cut = 0.15, say. As already observed in connection with the 
Xq F values, sets LO P and NLO P give the best description of the data. The contribution 
due to gluon fragmentation, which only enters at NLO, is insignificant, below 1%. The 
contribution due to the first four quark flavors is mostly concentrated at low x and is also 
very small. For x > x cnt , it makes up less than 1% of the total integrated cross section. 
It is interesting to study the b — > B branching fraction, 

B b (n) = I dxD b (x,fx), (4) 

where, for reasons explained above, we have introduced a lower cutoff at x cut = 0.15. 
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In Table |2], we present the values of B b (fi) at threshold /i = 2m?, and at the Z-boson 
resonance \i = M z for the various FF sets. As expected, B b (n) is rather stable under 
the evolution from 2m b to Mz- The values of Bb(M z ) are consistent with the input 
f{b — > B) — 0.405 ± 0.035 ± 0.045 |29| which was used to scale the experimental data 



points JTTj] so as to obtain the fully normalized cross section. 

Another quantity of interest is the mean B to b momentum fraction, 

1 /-i 

(x) b (ii) = dxxD b (x,fx). (5) 

-D&(/t) ■'Scut 

Table ^|also contains the values of (x) b (p,) at [A — 2m b and Mz evaluated with the various 
FF sets. The differences between sets S, P, and B on the one side and between LO and 
NLO on the other side are small. As fi runs from 2m b to Mz, (x)b(fi) decreases from 
approximately 0.8 down to below 0.7. This is a typical feature of the /i evolution, which 
generally softens the FF's. Our values of (x) b (Mz) can be compared with the experimental 
result reported by OPAL |TT| , 



(x)b(M z ) = 0.695 ± 0.006 ± 0.003 ± 0.007, (6) 

where the errors are statistical, systematic, and due to model dependence, respectively. 
Our results in Table |2| are in reasonable agreement with Eq. (^[). In connection to this, 
we remark that Eq. (|j) is not directly obtained from the measured distribution, which 
would be difficult to do, since there are no data points below x = 0.2. To extrapolate to 
the unmeasured region, OPAL uses four different models which describe the primordial 
fragmentation of b quarks inside their MC simulation. Equation is actually determined 
from the MC fits to the measured data points. Obviously, the quoted error for the model 
dependence can only account for the specific model dependence inside their particular 
MC approach, and need not be characteristic of the absolute model dependence. A 
rather model-independent fit to the x distribution, including a MC estimate for the region 



x < 0.2, leads to (x)b(M z ) = 0.72 ± 0.05 where the error is only statistical and 



does not account for the uncertainty due to the extrapolation. Our results in Table |] 
are somewhat smaller than this value and are barely consistent with the experimental 
error given above. Nevertheless, we believe that our results in Table |2| are in reasonable 
agreement with the independent determinations of (x)b(Mz) quoted in Ref. |11 . 



3 _B-meson production in pp collisions 

In this section, we compare our LO and NLO predictions for the cross section of inclusive 
B + /B° production in pp collisions at the Tevatron (y^ = 1.8 TeV) with recent data from 
the CDF collaboration ||. These data come as the px distribution da/dpr integrated 
over the central rapidity region |?7| < 1 for pt values between 7.4 and 20 GeV. They are 
normalized in such a way that they refer to the single channel pp — > B + + X. In the case 
of run 1A, where both B + and B° mesons were detected, the respective cross sections 
were averaged, i.e., their sum was divided by a factor of two. 
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Our formalism is very similar to Ref. ||30|| , where inclusive light-meson production in 
pp collisions was studied in the QCD-improved parton model. The relevant formulas and 
references may be found in Ref. and we refrain from repeating them here. We work 
at NLO in the MS scheme with rif — 5 massless flavors. For the proton and antiproton 
PDF's, we use set CTEQ4M {HJ with = 202 MeV. We evaluate a s from the two-loop 
formula with this value of A^. We recall that the evolution of the FF sets NLO S, NLO P, 

and NLO B is performed with A^ = 227 MeV, which is very close to the above value. 
We identify the factorization scales associated with the proton, antiproton, and the B 
meson and collectively denote them by Mf. We choose renormalization and factorization 



scales to be [i — Mf = 2iiit, where my = yp^ + mj, is the 5-meson transverse mass. 
Whenever we present LO results, they are consistently computed using set CTEQ4L |yj 
of the proton and antiproton PDF's, our LO sets of 5-meson FF's, the one-loop formula 



(5) 

for a s with A^ = 181 MeV |pl| , and the LO hard-scattering cross sections. We adopt 



the kinematic conditions from Ref. ||. Since we employ Di,(x,fi) both for the b and b 
quarks, the resulting cross section corresponds to the sum of the B + and B~ yields. Thus, 
it needs to multiplied by a factor of 1/2, in order to match the cross section quoted in 
Ref. 1. 

First, we consider the px distribution da/dpr integrated over the rapidity region \rj\ < 1 
as in the CDF analysis ||. In Fig. @(a), we compare the CDF data || with the LO and 
NLO predictions evaluated with our various sets of 5-meson FF's. The NLO distributions 
fall off slightly less strongly with increasing p? than the LO ones. The results for sets 
LO S, LO P, and LO B almost coincide. The same is true of the results for sets NLO S, 
NLO P, and NLO B. This means that the details of the b — > B fragmentation is tightly 
constrained by the LEP data, and that the considered variation in the functional form 
of the 6-quark FF at the starting scale has very little influence on the pt distribution. 
Henceforth, we shall only employ sets LO P and NLO P, which yielded the best fits to 
the OPAL data [IT]. We observe that our prediction agrees very well with the CDF data, 
within their errors. This is even true for the data point with smallest px, Pt — 7.4 GeV, 
where the massless approach is presumably not valid any more. It should be emphasized 
that the NLO prediction reproduces both the shape and the absolute normalization of the 
measured pt distribution, while the previous investigations mentioned in the Introduction 
fell short of the data by roughly a factor of two. 

The CDF collaboration has not yet presented results on the 77 distribution of the 
produced B mesons, which would allow for another meaningful test of the QCD-improved 
parton model endowed with 5-meson FF's. Anticipating that such a measurement will 
be done in the future, we show in Fig. |2|(b) the r] dependence of d 2 a/dt]dpT evaluated with 
sets LO P and NLO P at p T = 13.4, 17.2, 20, and 30 GeV. The first three of these p T 
values are among those for which CDF performed measurements of da/dpT [§]• Since the 
r\ spectrum is symmetric around rj = 0, we only consider 77 > in Fig. 0(b). As expected, 
the cross section falls off with 77 increasing from zero up to the kinematic limit, which 
depends on px- 

In order to assess the reliability of our predictions, at least to some extent, we now 
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investigate the scale dependence of the cross section considered in Fig. 0(a). To this end, 
we introduce the scale factor £ such that \i = Mf = 2^m T . In Fig. [| the £ dependence 
of da/dpT is displayed for pt = 13.4, 20, and 30 GeV. The calculation is performed with 
sets LO P and NLO P. For the two highest pr values, pr = 20 and 30 GeV, we observe 
the expected pattern. The LO results for da/dpr essentially decrease with £ increasing, 
whereas the NLO results are rather £ independent and exhibit points of horizontal tangent 
close to £ = 1. Furthermore, the LO and NLO curves intersect near these points. Thus, 
the scale choice £ = 1 is favoured both from the principles of minimal sensitivity 



and fastest apparent convergence |33[ . These observations reassure us of the perturbative 
stability and the theoretical soundness of our calculation in the upper pt range. For 
Pt = 13.4 GeV, the NLO prediction of da/dp T shows a stronger scale dependence, in 
particular, when the scale is drastically reduced. If we limit the scale variation to the 
interval 1/2 < £ < 2, which is frequently considered in the literature, the NLO cross 
section still varies by a factor of 1.56, to be compared with 1.15 at pt = 20 GeV. We 
hence conclude that, below px = 13.4 GeV, our NLO predictions should be taken with a 
grain of salt. The dents in the curves for p? = 13.4 GeV appear at the value of £ where 
Mf = fj, . This is because we identify D b (x,Mf) = D b (x,fj, ) if Mf < fi , i.e., the FF's 
are frozen below their threshold. 

We must also remember that, for pt values comparable to m b , the massless-quark 
approximation ceases to be valid, since terms of order mf/p T are then not negligible 
anymore. For pt = 13.4 GeV and 20 GeV, we have m\jp\ = 0.14 and 0.063, respectively, 
so that the massless approximation should certainly be valid for p? = 20 GeV. On the 
other hand, for p T = 20 GeV, we have a s ln(pf,/m|) = 0.42, assuming that a s = 0.15, so 
that the NLO calculation in the massive scheme, where these logarithmic terms are not 
resummed, should then already be inadequate. From these considerations, we conclude 
that our predictions should be fairly reliable for p T ^ 15 GeV. 



4 Conclusions 

In this paper, we considered the inclusive production of single B mesons in the QCD- 
improved parton model endowed with nonperturbative FF's. We chose to work at NLO 
in the pure MS factorization scheme with five massless quark flavors. This theoretical 
framework is known to lead to an excellent description of a wealth of experimental in- 
formation on inclusive light-hadron production in different types of high-energy reactions 



27| . |34]| . It is thus expected to also work well in the case of B mesons provided that the 
characteristic mass scale M of the process by which they are produced is large compared 
to the 6-quark mass. Then, the large logarithms of the type a s \n(M 2 /ml) which are 
bound to arise in any scheme where bottom is treated as a massive flavor get properly 
resummed by the AP evolution, while the omission of the terms suppressed by powers 
of m\jM 2 is a useful approximation. The criterion M 3> m;, is certainly satisfied for 
e + e~ annihilation on the Z-boson resonance, and for hadroproduction of B mesons with 
Pt ^> irib- Owing to the factorization theorem, the FF's are universal in the sense that 
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they just depend on the produced hadrons and the partons from which they sprang, but 
not on the processes by which the latter were produced. Thus, the theoretical framework 
adopted here is particularly suited for a consistent description of LEP1 and high-p^ Teva- 
tron data of inclusive 5-meson production. By the same token, a massive calculation at 
fixed order would be inappropriate for this purpose. 

Our procedure was as follows. We determined LO and NLO 5-meson FF's by fitting 
the fractional energy distribution of the 5-meson sample collected by the OPAL collabora- 



tion at LEP1 |TT|. In order to get some handle on the theoretical uncertainty, we adopted 
three different functional forms for the b — > B FF at the starting scale, which we took to 
be /j,o = 2m& = 10 GeV. The ansatz proposed by Peterson et al. yielded the best LO 
and NLO fits, with Xdf = 0-67 and 0.27, respectively. The e parameter, which measures 
the smearing of the Peterson distribution, came out as 0.0126 and 0.0198, respectively, 
i.e., more than twice as large as the standard value e = 0.006 determined by Chrin [ 10f 



more than a decade ago, before the LEP1 era. In this connection, we should emphasize 
that the results for the fit parameters, including the value of e, are highly scheme depen- 
dent at NLO, and must not be naively compared disregarding the theoretical framework 
to which they refer. The b — ► B branching fraction and the mean B to b momentum 
fraction evaluated from the resulting FF's after the evolution to [i = Mz turned out to be 
in reasonable agreement with the model-dependent determinations by OPAL |IT |. Using 



our FF's, we made theoretical predictions for the inclusive hadroproduction of single B 
mesons with large px- We found good agreement, both in shape and normalization, with 
the pt distribution recently measured in the central rapidity region by the CDF collab- 
oration at Fermilab ||. From the study of the scale dependence of the LO and NLO 
calculations, we concluded that our results should be reliable for p T ^ 15 GeV. To our 
surprise, the central prediction, with scales /i = Mf = 2rriT, also nicely agreed with the 
CDF data in the low-pr range, where the massless scheme is expected to break down. We 
recall that the massive NLO calculation with traditional Peterson fragmentation [l(| was 
found to fall short of these data by a factor of two. It would be interesting to also test 
the predicted rj distribution against experimental data. 
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TABLE CAPTIONS 



Table 1: Fit parameters for the b —> B FF's according to sets S, P, and B at LO and NLO 
and respective values of y 2 per degree of freedom. All other FF's are taken to be zero at 
the starting scale /io = Ira^ = 10 GeV. 

Table 2: b — > B branching fractions and mean B to b momentum fractions evaluated from 
Eqs. (H) and (f|), respectively, at the starting scale and at the Z-boson resonance using 
the various FF sets. 



FIGURE CAPTIONS 



Figure 1: The cross sections of inclusive B + / B°-meson production in e + e~ annihilation 
at ^Js = Mz evaluated with sets (a) LO S and NLO S, (b) LO P and NLO P, and (c) 
LO B and NLO B are compared with the OPAL data fll| . 

Figure 2: (a) The cross section da/dpx of inclusive i? + /i? -meson production in pp colli- 
sions with -v/i =1.8 TeV, integrated over \t]\ < 1, is compared with the CDF data f|. The 
predictions are calculated at LO and NLO with sets S, P, and B. (b) The cross section 
da I drjdpT at fixed values of pt evaluated with sets LO P and NLO P. 

Figure 3: Scale dependence of the cross section da/dpr, integrated over \r)\ < 1, at fixed 
values of pt- The predictions are calculated at LO and NLO with set P. 
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set 


N 


a 


P 


e 


r 


Xdf 


LO S 


56.4 


8.39 


1.16 






0.80 


NLO S 


79.4 


8.06 


1.45 






1.21 


LO P 


0.0952 






0.0126 




0.67 


NLO P 


0.116 






0.0198 




0.27 


LO B 


0.308 








0.121 


2.50 


NLO B 


0.280 








0.156 


1.66 



Table 1 



set 


B b (2m b ) 


B b (M z ) 


(x) b (2m b ) 


(x) b (M z ) 


LO S 


0.425 


0.411 


0.813 


0.697 


NLO S 


0.384 


0.370 


0.787 


0.672 


LO P 


0.448 


0.431 


0.787 


0.677 


NLO P 


0.405 


0.388 


0.758 


0.650 


LO B 


0.460 


0.442 


0.768 


0.663 


NLO B 


0.416 


0.398 


0.739 


0.635 



Table 2 
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